Staggered Flux Phase in a Model of Strongly Correlated Electrons 



J. B. Marston and J. 0. Fjaerestad 
Department of Physics, Brown University, Providence, RI 02912-1843 



A. Sudb0 

Department of Physics, Norwegian University of Science and Technology, N-7491 Trondheim, Norway 



(N 
O 
O 
(N 

C 



h 

c3 



i 

c 

o 
o 



> 

00 

oo 

o 

O 

c3 



c 

o 
o 



X 
S3 



(Dated: February 11, 2002) 

We present numerical evidence for the existence of a staggered flux (SF) phase in the half- 
filled two-leg t-U-V-J ladder, with true long-range order in the counter-circulating currents. The 
density-matrix renormalization-group (DMRG) / finite-size scaling approach, generalized to de- 
scribe complex-valued Hamiltonians and wavefunctions, is employed. The SF phase exhibits robust 
currents at intermediate values of the interaction strength. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.30.+h, 74.20.Mn 



The zero-temperature phase diagram of the two- 
dimensional Hubbard model and its various extensions 
remains poorly understood. An intriguing possibility 
that has been the focus of considerable attention is 
known variously as the "orbital antifcrromagnet" Jl| p| , |) , 
the "staggered flux" (SF) phase @, |, §, % |, ff! or 
the "d-density wave" 0, The state breaks both 

time-reversal and lattice translation symmetries; another 
phase, the "circulating current phase" |L2| , is similar 
but preserves translational symmetry. Such "hidden" 
forms of order could arise in the pseudogap phase of 
the cuprate superconductors, implying the existence of 
a quantum critical point and possibly also non-Fermi liq- 
uid behavior jllj The SF phase competes against 
other better-known phases such as charge-density waves 
(CDW), spin-density waves (SDW), superconductivity 
(BCS), stripes, and Fermi liquids. It is important to 
ascertain whether or not such order can really occur in 
sensible microscopic models of correlated electrons. In 
this Letter we study numerically the simplest possible 
system that could support such a phase, namely an ex- 
tended version of the Hubbard model on a two-leg ladder 
(see Fig. ^). Our approach extends to arbitrary inter- 
action strength earlier analytical work by two of us that 
showed that the SF phase arises at weak-coupling |T^| . 

Ladder systems are of interest both because they 
are easier to understand theoretically than full two- 
dimensional lattices, and also because they are realized in 
nature |l4|]. Although ladders are one-dimensional, true 
long-range order (LRO) in the orbital currents is possible 
at zero temperature because the currents break discrete 
symmetries. This is in contrast to spin-density wave or 
BCS superconducting order which break continuous sym- 
metries and which therefore can at most exhibit quasi- 
LRO with power-law decay. 

Away from half-filling, numerical studies of the t-J 



model |l6| and various extensions of it|l7| on the two- 
leg ladder found no evidence for a SF phase. However, for 
weak interactions, tendencies towards SF ordering (not 
true LRO, but rather power-law decay of current-current 
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FIG. 1: Upper diagram: A two-leg ladder of length L = 8. 
The left and right blocks each retain Hilbert spaces of dimen- 
sion M . To implement the DMRG / finite-size algorithm, the 
two sites belonging to each rung are paired together into a sin- 
gle site with a Hilbert space of complex dimension D = 16. 
Lower diagram: Source current jsF applied to the left-most 
rung induces currents in the interior of the ladder (dashed 
circles). A source of Cooper pairs {jbcs), a local chemical 
potential (jcDW/Sp), an d a local magnetic field (jsmv) are 
added at the ladder ends to induce BCS, SP, CDW, and SDW 
order in the central region of the ladder, which is then moni- 
tored as the ladder length increases. 



correlation functions) have been found analytically both 
for the spinless|I^| and spinfuljl^] cases (see also Ref. pp[ ). 

Here we will consider the half-filled two-leg ladder. We 
study an extended Hubbard model, the t-U-V-J model, 
defined on the two-leg ladder as follows: 
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(1) 



Here cj^ creates an electron of spin a =t, J. on site i of 
leg A = t, b of the ladder. Operators n,,A = cj^c^A.c 

and Si y \ = |c| " (?f c^a,/3 are respectively the number and 
spin at site (i, A). There is an implicit sum over repeated 
raised and lowered spin indices. Particle-hole symmetry 
implies that the chemical potential fi = when the sys- 
tem is half-filled, with (n^x) = 1 on each site. 

For Vj| = Jy = and Jj_ = 4(f7 + VI) the model 
has SO(5) symmetry The phase diagram of the 
SO(5)-symmetric model with t±_ = in was investigated 
in Ref. ^2| in the weak-coupling limit using a perturba- 
tive renormalization-group (RG) analysis combined with 
bosonization. One of the phases of the SO (5) model, lo- 
cated in the U — V± plane between the lines V± = —2U 
and V± w —5.7U with V±_ > and U < 0, was identi- 
fied as having spin-Peierls (SP) order. However, a recent 
bosonization study by two of us [l3| revealed that this 
phase in actuality exhibits no dimerization, but instead 
is the SF phase with true LRO in the orbital currents. 
For a correct understanding of the problem it was essen- 
tial to treat carefully the "Klein factors" that must be 
introduced to maintain anti-commuting statistics of the 
(bosonized) fermionic degrees of freedom. 

Previous RG studies have shown that in weak cou- 
pling a rather generic ladder model flows to a manifold 
with SO (5) symmetry for a range of values of the model 
parameters p2], p3[ . In agreement with these results, we 
have shown that there are RG flows towards the SF phase 
also when not too large non-SO(5)-invariant terms (such 
as J|| and Vj|) are added to the SO(5)-modelp^]. 

Weak-coupling RG is unreliable at intermediate in- 
teraction strengths. Instead we use the "infinite- 
size" version of the density-matrix renormalization-group 



(DMRG) algorithm^! to search for the ordered currents. 
The half-filled ladder is expected to be fully gapped over 
a large portion of the phase diagram; consequently (as 
shown below) the infinite-size algorithm is sufficiently ac- 
curate to obtain well-converged results. Each site of the 
ladder has a Hilbert space of dimension 4, as the site 
can either be unoccupied, have a single electron of ei- 
ther up or down spin, or be doubly occupied. To employ 
the DMRG algorithm, we group pairs of sites connected 
by a rung into a single site of Hilbert space dimension 
D = 4 x 4 = 16 (see Fig. [l]). Errors in the calculation 
of observables introduced by the DMRG truncation of 
the Hilbert space can be systematically reduced by in- 
creasing the size M of the Hilbert space retained in each 
of the two outer blocks up to limits imposed by com- 
puter memory and speed. The largest Hilbert space we 
use has M = 150 with corresponding complex dimension 
M x D x D x M = 5, 760, 000. 

Questions of the spontaneous formation of order are 
addressed by the combined DMRG / finite-size scaling 
approach described in some detail in Ref. [2(| Quantum- 
critical points have been studied with the method, and 
critical exponents have been obtained at percent-level 
accuracy p6|, |7j. In the present case we induce sym- 
metry breaking by applying a source current jsf to the 
left end rung of the ladder. Terms that induce Cooper- 
pair formation by breaking U(l) particle-number sym- 
metry {jbcs), induce CDW and SP order by breaking 
lattice reflection symmetries (jcDW/Sp): an d induce a 
SDW through the application of a local magnetic field 
Usdw), are also added to the right end of the ladder 
(see Fig. 0): 



H -> H + j SF * it ± (c\° t cx,b,o - H.c.) + jbcs * (4!t c i,6 " c % c L,b + H - c -) 

+ JCDW/SP * n L,t + jsDW * SI . (2) 



r 



In addition to these explicit source terms, we note that 
the open boundary conditions on the ladder also act as 
0(1) source terms for a columnar dimer pattern|2l|. As 



the Hamiltonian is complex- valued when jgp ^ 0, we 
generalize the DMRG code to handle complex-valued 
wavefunctions and reduced density matrices. This gener- 
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alization comes at the cost of doubling the required mem- 
ory, and slowing down both the sparse and the dense ma- 
trix diagonalization, but the cost is offset by the advan- 
tage that now we can access fully ordered ground states, 
well beyond linear response theory. 

We calculate the expectation value of the current oper- 
ators 2t± lm{c\ a t Ci^,a} and 2£y Imjcj^ x Ci^ a } respec- 
tively on the central rungs and links of the ladder, check- 
ing that current is conserved (Kirchhoff 's junction rule) 
at the vertices in the central region; deviations are typi- 
cally of order 10~ n in units where t± = t\\ = 1. We also 
monitor the Cooper pair amplitude, (c|jc|^— c} \c\ \); the 
average magnitude of deviations in the electron hopping 
amplitudes {c\\cj_ v,<r + H.c.) from their mean value; de- 
viations in the electron occupancy from one, (n^\ — 1); 
and the local spin density, (S* x ). 

First consider the pure minimal Hubbard t-U model 
with t± = i|| = 1, U = 4, and Vj_ = V\\ = J± = J\\ = 0. 
As shown in the semilog plot of Fig. ||(a) , applied source 
currents of jsF — 1 and 0.01 induce currents on the four 
central rungs of the ladder that decrease exponentially as 
the ladder is enlarged via the DMRG algorithm, in agree- 
ment with the weak-coupling RG calculation (22). Like- 
wise we find no instabilities towards BCS, SP, CDW or 
SDW order. Instead the phase is a fully-gapped insula- 
tor qualitatively the same as that found in the Heisenberg 
antifcrromagnet . 

Next consider the effect of turning on interactions J± 
and V± along the rungs of the ladder. According to 
the phase diagram in Ref. ^2|, for the case of exact 
SO(5) symmetry with J± = 4(U + V±) the SF phase 
should occur p| in the weak-coupling limit for U < 
at V±/U = -9/4. Fig. | shows that the SF phase 
does indeed arise at these ratios for intermediate inter- 
action strength: For U = —0.4, V± = 0.9, and J± = 2, a 
small source current induces orbital currents in the cen- 
tral region of the ladder that grow with increasing lad- 
der length. The currents saturate to the same non-zero 
asymptotic value in the limit of large ladder length, re- 
gardless of the size of the source terms. There is good 
convergence with increasing block Hilbert space dimen- 
sion M] the asymptotic value of the rung current for the 
most accurate truncation of M — 150 is 0.8047 and dif- 
fers by only 0.1% from the value obtained with M = 100. 
The currents alternate in sign from rung to rung as ex- 
pected in the SF phase. There are no other instabil- 
ities: As shown in Fig. | the BCS, SP, CDW, and 
SDW order parameters all vanish exponentially as the 
ladder is enlarged via the DMRG algorithm, consistent 
with the weak-coupling prediction. The excitation spec- 
trum remains fully gapped; consequently the SF phase 
exists over a range of parameter space well beyond the 
region of exact SO (5) symmetry. For example, robust 
currents arise for U — —0.5, V± = 0.75, and J± = 2, 
and also when interactions along the links are turned on, 
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FIG. 2: Semilog (a) and linear (b) plots of the magnitude of 
the induced current averaged over the four central rungs of the 
ladder. Dotted lines are for calculations with M = 50, dashed 
lines are for M — 100, and solid lines are for M = 150. Results 
are shown for the case of the SO(5) ladder with t± = t|| = 1, 
U = -0.4, Vx = 0.9, J± = 2, and V n = J\\ = 0. Also 
shown on the semilog plot is the case of the minimal Hubbard 
model with (7 = 4 and jsF = 0.01 (lower line that drops 
exponentially) and j S F = jscs = Jcdw/sp = jsmv = 1 
(upper lines that drop exponentially) . In both plots, the lower 
set of SO (5) curves is for a small source current jsF = 0.01 
applied to the ladder's left edge. The upper set is for large 
source currents jsF = Jbcs = Jcdw/sp = JSdw = I- Note 
the convergence to the same asymptotic value of current as 
the thermodynamic limit of long ladder length is approached. 
The asymptotic current for M — 100 differs by 0.6% from 
that at M = 50; the value for M = 150 differs by only 0.1% 
from that at M — 100 demonstrating good convergence. 



U = -0.5, V x = Vj| = 0.75, Jx = 2, and Jj| = 0.5. 

In the large- Jj_ limit, with ratios U/Jx = —1/4 and 
Vx/Jx = 3/8 kept fixed, the ground state of the ladder 
consists of a direct product of spin-singlet dimers on each 
of the rungs, with no broken symmetries pl| . We have 
verified that the orbital currents cease when the interac- 
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FIG. 3: Semilog plot of various order parameters in the SF 
phase (see text). Edge source terms jbcs = jcDw/sp = 
jsDW — jsF ~ 1 induce modulations in corresponding ob- 
servables at the center of the ladder. Dotted, dashed, and 
solid lines correspond to M — 50, 100 and 150 respectively. 
Plotted in the figure, from highest to lowest in magnitude, 
are the order parameters for CDW, SP, BCS, and SDW ten- 
dencies. See text for the definition of these order parameters. 
Each of the expectation values decays exponentially with in- 
creasing ladder length. (The SP order parameter eventually 
levels off, but this residual order vanishes as the block size 
AI is increased, showing that it is an artifact of the Hilbert 
space truncation.) Similar results are obtained upon applying 
either large or small source terms, one at a time. Thus there 
is no tendency towards BCS, SP, CDW or SDW order in the 
SF phase. 

tion strength is increased further to J± > 2.5. We are 
carrying out a more extensive investigation of the phase 
diagram of model Eq. [I] by using the weak-coupling RG 
equations in combination with the DMRG / finite-size 
scaling method |24j]. 

It would be interesting to study the effects of doping 
away from half-filling. Numerical study of the doped SF 
phase would require the use of the more accurate "finite- 
size" DMRG algorithm, as either gapless excitations, or 
rich spatial structures such as stripes, are expected to 
occur pgj]. One question to be answered is whether or 
not LRO in the currents can still arise at commensurate 
hole concentrations. 

In summary we have shown that a phase of strongly 
correlated electrons exists in which currents form sponta- 
neously in the half-filled two-leg ladder. At half-filling the 
SF phase is fully gapped, and thus is an insulator of the 
Mott-Hubbard type. Our results, suitably generalized to 
two or three spatial dimensions, may have application 
to several stoichiometric compounds that were recently 
proposed to be in the SF phase |^9|. It is also intriguing 
that the SF phase lies in between a checkerboard CDW 
phase, and the D-Mott phase jl3|, |2^, as charge seg- 



ration into stripes and d-wave superconductivity are two 
phases that occur in the high-T c cuprates. 

Note added: U. Schollwock has recently used the 
"finite-size" version of the DMRG algorithm, generalized 
to complex- valued wavefunctions, to study the SO (5) lad- 
der with the same parameters as Fig. |2| The value of the 
saturated rung current that he finds agrees quantitatively 
with our result. 
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